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1. Introduction 

In a recent paper (|astro-ph/0207561] , posted on July 25, 2002), Fromang & Balbus 
(2002, hereafter FB02) claim that the minimum energy configuration of a rotating fiuid 
depends critically on whether or not any part of the fiuid rotates supersonically. They 
O ' state that "a uniformly rotating barotropic fiuid in an external potential attains a true 
^ ■ energy minimum if and only if the rotation profile is everywhere subsonic. If regions of 
supersonic rotation are present, fiuid variations exist that could take the system to states 
of lower energy." 
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We present here a simple counter-example which demonstrates that the stablity cri- 
terion presented in FB02 is incorrect. 



^ ! 2. Minimum Energy Equilibria 

^ I We first repeat part of their analysis. Following FB02, we consider a barotropic fiuid 

OO ' in which the pressure is a unique function of the density. To further simplify matters, and 

^ ■ to be consistent with the analysis in the next section, we consider a polytropic equation of 

. state in which the pressure P and the density p are related as 

p.: p = Kp\ (1) 

o , 

^ , Here iiT is a constant throughout the fiuid and is not allowed to vary in any perturba- 
^ \ tion. Thus we exclude thermal effects. Furthermore, in line with FB02, we consider only 
perturbations that are adiabatic, which excludes convective instabilties. 

For such a fiuid, the internal energy per unit volume e, the enthalpy per unit mass H, 
and the adiabatic sound speed a, are given by 



P H=^P, a^=2f. (2) 
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We assume that the mass M and the total angular momentum J of the system are 
held fixed: 

M = J pdV = constant, (3) 

J = y pRvdV = constant, (4) 

where R is the cylindrical radius as measured perpendicular to the rotation axis and v{R) 
is the azimuthal velocity of the gas at radius R. Note that, as do FB02, we take account 
throughout only of the azimuthal component of the velocity. The implication of this is 
that we restrict stability considerations to axi-symmetric perturbations alone. 
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We are interested in minimizing the energy, E, of the system, 



E = I (^^pv^ + e + pcjj dV, (5) 

where (f) is the (time-independent) external potential in which the fluid is confined. (Fol- 
lowing FB02, we neglect self-gravity of the fluid.) The minimiziation of E subject to the 
constraints (3) and (4) is a standard problem that can be solved by variational methods. 

I = E -CM -QJ = J (^^pv'^ + 6 + p(l)-Cp-^pRv^dV, (6) 

where C and il are Lagrange multipliers. To find an extremum of the energy, we set the 
first order variation of / with respect to changes in p{R) and v{R) to zero: 



5/ = y Qw^ - QRv + H + <t)-C^Sp + p{v- nR)dv 



dV = 0, (7) 



where we have used equation (1) to set de/dp — H. Since SI should vanish for arbitrary 
choices of Sp{R) and 5v{R), the coefficients of 5p{R) and 5v{R) in the integrand of (7) 
should individually vanish at each R. This gives the following two conditions: 

V = nR, (8) 

-^fi^i?2 + H + (f) = C^ constant. (9) 

The first condition states that a polytropic fluid in an energy-extremum configuration must 
be in a state of uniform rotation, a well-known result (e.g. Landau & Lifshitz). The second 
condition is nothing other than the condition of hydrostatic balance, as can be verified by 
taking the gradient of equation (9). We thus find that a polytropic fluid with fixed M and 
J and with an extremum in E has uniform rotation and is in hydrostatic equilibrium. 

The analysis is identical to that in FB02, with the sole difference that we have spe- 
cialized to the case of a polytropic equation of state. Following the above steps, FB02 
proceed to take second differences of the energy, on the basis of which they reach the con- 
clusions stated above in §1. We do not repeat their analysis since their final conclusions 
are incorrect, as we now demonstrate by means of a counter-example. 



3. A Simple Counter-Example 

We consider a polytropic fluid rotating around the z axis in a quadratic potential 

(P{R,z)^^nl{R^ + z^). (10) 

We assume that the fluid is in an energy extremum state which, by the analysis of the 
previous section, means that the fluid rotates with a constant angular velocity Q and is 
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in hydrostatic equilibrium, satisfying equation (9). Let us assume that the fluid extends 
from i? = to a maximum radius Rq on the equatorial plane. (Such a maximum radius 
always exists for 7 > 1.) At the point {R,z) = (Rq^O), the enthalpy of the gas vanishes. 
This enables us to determine the constant C in equation (9), using which we can show that 



H{R,z)^^{Qj,-n') 



9? 



(11) 



Not surprisingly, the fluid takes the form of a flattened spheroid, with equatorial radius 
Rq and "vertical" radius 



i?o. (12) 



From the expression for we may calculate the density p of the fluid, 

^ r(7-i)^r 
P(i?,^)=[^^^_ , 

where we have deflned the polytropic index n in the usual way: 



(13) 



n=^—. (14) 
7-1 ^ ' 

The vertically integrated surface density of the fluid is then 

S(i?) = C{Rl - i?2)«+V2^ (15) 

_ r(i/2)r(n+ 1) / 7 - l y (^1. - ^'r^''' .^g. 

r(n + 3/2) ) Qk ' ^ ^ 

Through equation (16), the angular velocity f2 deflnes the parameter C. Thus this 
simple model has two free parameters, Rq and Q. For a given value of Q we may determine 
the radius Rq by applying the mass constraint 

M = J E{R)27TRdR = ^^^^i^f Rl""^^- (17) 

We then determine VL by applying the angular momentum constraint, 

J = I J:{R)QR'^2nRdR = -. ^—-MnRl (18) 

J (n + 5/2) 

From equations (16) and (17) we see that for a given mass, Rq increases with increasing 
O, and thus from equation (18) it is clear that J is also a monotonically increasing function 
of ft. We conclude that for given values of M and J, there is one and only one energy 
extremum conflguration. 
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All possible fluid configurations in the model have a positive energy, since all the 

terms in equation (5) are positive. Therefore, we are guaranteed that the system has a 
well-defined global energy minimum. Since we have shown that for given values of mass 
M and angular monentum J, there is only one configuration that is an energy extremum, 
that particular configuration must correspond to the absolute energy minimum state and 
must be stable. 

In this simple analytical model, the orbital velocity v{R) increases monotonically with 
radius i?, from t; = 0ati? = 0toa maximum aX R = Rq, while the sound speed a{R, z) 
decreases from a maximum at -R = 0, ^ = to a = all over the surface (and specifically 
at the outer edge R = Rq^z = Qi): 

v{R) = nR, a\R, ^) - (7 - l)H{R, z)^[l-^^-^. (19) 

It is evident that the model includes both subsonic and supersonic zones, and should 
therefore be unstable according to the criterion given by FB02. And yet, as proved above, 
the fiuid is in an energy minimum state and is stable. The model is thus an explicit 
counter-example to the result claimed in FB02 — supersonic rotation does not imply that 
the fiuid has neighboring states with a lower energy. 

4. A Final Comment 

In astrophysics, one is generally interested in fiuids with free surfaces at zero pressure 
rather than fiuids confined between walls. Since FB02 claim that the only stable barotropic 
configurations are those that rotate uniformly and are subsonic at all radii, it is interesting 
to ask whether such fiuid configurations are at all possible in the first place. The answer 
is that, if we leave out the strictly non-rotating case (fi = 0) and if we restrict ourselves to 
fiuids with polytropic indices n > 0, there are no isolated configurations that meet FB02's 
requirements. (We thank Martin Rees for pointing this out.) 

For instance, if the fluid extends to an infinite radius (as can happen for an isothermal 
gas), then the rotation velocity v{R) — > oo at large R whereas the sound speed remains 
finite. The fiuid at large radius is thus guaranteed to be supersonic. On the other hand, 
if the fiuid has an edge at a finite radius, as in the simple model described in §3, then v 
is finite at the edge whereas the sound speed goes to zero. Once again, the fiuid near the 
outer edge rotates supersonically. 

Thus, all uniformly rotating configurations for a wide class of fluids have supersonic 
motions. Since, surely, all of these fluids must possess stable minimum energy equilibria, 
the result claimed by FB02 must be incorrect. The example discussed in §3 is a speciflc 
example that explicitly demonstrates this. 
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